Cosmic Sparks from Superconducting Strings 
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We investigate cosmic sparks from cusps on superconducting cosmic strings in light of the recently 
discovered millisecond radio burst by Lorimer et al\j^. We find that the observed duration, fluence, 
spectrum, and event rate can be reasonably explained by Grand Unification scale superconducting 
cosmic strings that carry currents ~ 10 GeV. The superconducting string model predicts an event 
rate that falls off only as S^ 1 ^ 2 , where S is the energy flux, and hence predicts a population of very 
bright bursts. Other surveys, with different observational parameters, are shown to impose tight 
constraints on the superconducting string model. 



The discovery of a radio burst ( "spark" ) was recently 
reported by Lorimer et al [l| in the Parkes survey and 
analysis of the dispersion measure suggests that the 
source is of cosmological origin. A thorough examina- 
tion of the observation has supported this conclusion Q • 
Conventional astrophysical sources are not known and 
neither has a host galaxy for the event been identified. 

If more cosmic sparks are observed and found to be 
extra-galactic, it would indicate an exotic cosmological 
process. Superconducting cosmic strings are a possible 
exotic source of electromagnetic phenomena in the uni- 
verse and arise quite naturally in particle physics mod- 
els 0], though less so in string theory 0]. Earlier work 
on cosmological signatures of superconducting cosmic 
strings has primarily foucussed on high energy emission 
in the form of particles and/or gamma rays [3J, [5|, l(| or 
synchrotron emission [3j . In this paper we argue that ra- 
dio emission may be a good way to look for superconduct- 
ing strings as they can cause observable sparks similar to 
the one seen by Lorimer et al. Even if further observa- 
tions discover astrophysical sources for observed sparks, 
the prediction that superconducting cosmic strings pro- 
duce radio sparks that are potentially observable is im- 
portant from the particle physics viewpoint, since their 
detection or absence may be used to constrain various 
fundamental models. 

Superconducting cosmic strings can be viewed as elas- 
tic, current-carrying wires, distributed in the cosmos as 
closed loops and infinitely long Brownian curves. The 
mass per unit length of a string will be denoted by fi = rj 2 
and the current by iq. Current-carrying strings oscillate 
under their own tension and radiate electromagnetically 
(e.g. see Q). The radiation is very strong from events 
such as "cusps" which are points on an idealized (zero 
thickness, no current) string that reach the speed of light 
for a brief instant [T3, EH ■ In a more realistic setting, 
the cusp gets cut-off due to the thickness of the string 
and due to the backreaction of the current and radia- 
tion. Nonetheless the radiation is very strong from local- 
ized regions ( "quasi-cusps" ) even in the realistic string 
case. The scenario we envision is that a curved section of 
string (or loop) of length L develops a cusp and beams 
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FIG. 1: String segment with current io at 3 different times, 
with a cusp at O, with velocity along z. The size of the curved 
section of string is L. The observer is at an angle Oo from the 
z-axis at a distance d. 



electromagnetic radiation in direction z. The observer is 
located at some large distance, d, from the location of the 
cusp and slightly off the z axis, at an angle 9q (Fig. [I]). 

The energy emitted from a cusp at angular frequency 
lu = 2ttv, at angle 9 from the beam direction, and ob- 
served at distance d (i.e. the fluence) is given by 



F = 



1 d 2 E UJ 
d 2 dtodQ 



w 2L 



-aujL8 3 



if 



jL9 3 > 1 (1) 



where a and b are constants that depends on the shape 
of the cusp and we will take a ~ 1, b ~ 1. (We also 
work in natural units throughout, so K = 1 = c.) This 
result is most easily obtained by expanding the result in 
Ref. [l(| for large u>L9 3 . If, instead, we have u>L9 3 < 1, 
the emitted energy falls off oc uj~ 2 I 3 . 

The total energy emitted from a loop in electromag- 
netic and gravitational radiation in one time period is 
-Etotaijoop ~ (nei r] + TGfi 2 )L = ^[iL where e RJ 0.3 
is the charge of the current carriers, G Newton's grav- 
itational constant, and n ~ 30, T ~ 100 are numerical 
factor determined by averaging numerically derived radi- 
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ation rates from a variety of loops. 

The temporal width of the emission seen from the cusp 
is given by the duration for which the observer is in the 
beam. The emission from the cusp is beamed but the 
beam is changing direction at a rate 9 ~ 1/L. So the 
beam changes its direction by 9q in a time L9q- Since, 
for large angles from the beam, the emission gets cutoff 
exponentially fast (Eq. ([1])), we have 9q ~ (wl) -1 ' 3 , and 



At ~ uj-^L 2 ' 3 



(2) 



The event rate depends on the number density of 
loops and is derived from numerical simulations of non- 
superconducting strings [g|. Note that the results for 
non-superconducting strings should also apply in the case 
that the current on the string is much smaller than the 
symmetry breaking scale 77. The simulations ignore radia- 
tive effects and find that the number density of loops of 
length between In and lo + dlo is dni ~ Adln / \%t 2 ) where 
A ~ 10 2 [HE Hi] (also see [H]). Due to radiation, the 
length I of a loop decreases with time, l(t) = l — j(t — ti), 
where ti is the time when the loop was born. So the num- 
ber density of loops of size / at time t is 



dni(t) 



Adl 



{I + -/t) 2 t 2 



(3) 



where the radiative shrinking has been included, and we 
have assumed that the loops were all born at some very 
early time so that t ^> £j. 

The spark observed by Lorimer et al was in the fre- 
quency interval (1.28,1.52) GHz and the central fre- 
quency is i>q — 1.4 GHz. Based on the dispersion mea- 
sure of the event, the observed event is constrained to lie 
within rcdshift of 0.1 — 0.3 and, for our estimates, we will 
assume that the event was located at z — 0.3, or at a 
comoving distance ~ 1 Gpc. The observed energy flux is 

ergs 



Sob 



30 Jy = 3 x 10" 



s-Hz 



(4) 



and the duration of the event has an upper bound At < 
Ato = 5 ms 

The observed pulse width cannot be used to estimate 
the intrinsic duration of the event because of scattering 
by the turbulent inter-galactic medium (IGM) The 
oberved time width, At oos , is a sum in quadratures of 
the intrinsic time width modified by cosmological redshift 
and the width due to scattering in the IGM: Ai b s = 
[(Ar obs ) 2 + (Ai cmit ) 2 ] 1/2 with 

At cmit = (1 + zy£gl*l* = co-^d + zf'H 2 ' 3 (5) 



The scattering time width at frequency v from an event 
at redshift z depends on the location of the scattering 
centers. For the case of a scattering screen close to the 
source, the width is @, Hij 



At ods (v,z) = At 



1 + z 
1 + zo 



/3+1 



1'0 



(6) 



where Ato = 5 ms, vq = 1.4 GHz, zq — 0.3, and 
(3 = —4.8. For comparison with the observed spark in 
the Parkes survey, we will use v = vq, but if comparing 
to other surveys it will be necessary to insert the appro- 
priate observational frequency. 

The observed power law fall off oc v~ A in the observa- 
tional frequency band can be fit by an exponential, and 
since the fluence is S oas At b s , 
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(7) 



To get an idea of the parameters needed of the su- 
perconducting string model, we first fit the observed 
spectrum ignoring redshift factors (which are small since 
z < 0.3). Using Eq. (Q} we find 



dkiF 
dlnuj 



;L9 3 = -4 



(8) 



The observation does not constrain the intrinsic du- 
ration of the event. However, if the intrinsic duration 
was 1 ms, we can equate the duration of the cusp event 
(Eqs. ([2])) to the observed duration to get lj^^L 2 / 3 = 
Ato = 1 ms, where ljq = 2-kvq ~ 10 10 s _1 . This gives 
L = Lul /2 (At) 3 / 2 ~ 3 s = 10 11 cm. If the intrinsic dura- 
tion were smaller than 1 ms, the length would be even 
smaller. This illustrates that we are considering signa- 
tures of loops that are very small on a cosmological scale. 

The energy flux from the cusp event is found from 
Eqs. (H|) and 1(5)). Equating to the observed value, Eq. |7)), 
gives 



i a ~ 10° 



10 n cm 
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GeV 



(9) 



where we have used 1 ergs/ (cm 2 — Hz) = 1 GeV 2 . Note 
also that the dynamics of the string will be dominated 
by the tension if 77 ^> 10 6 GeV. Hence it is a valid ap- 
proximation to ignore the current when discussing the 
dynamics of the string network. 

To estimate the rate of observed-sparkcr-like events, 
we use Eq. Then the event rate, denoted N, due to 
cusps on loops of size between I and l + dl, beamed within 
a solid angle dQ from us is 



dN - y 



Adl 



(l + "ft ) 2 t 2 4tt 



d ^dV 



(10) 



where the factors account for f c cusps per loop oscil- 
lation, the number density of loops, the beaming angle 
constraint, and the spatial volume. For Grand Unifica- 
tion scale strings, 77 ~ 10 14 GeV, we find 7/J0 ~ 10 9 s 3> 
I ~ 1 s, and we can rewrite (| 10[) as 



dN 



6 x — x sin 



x D dD 



(11) 
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where D = HqD c is the comoving distance to the loop 
in Hubble units, and C = 2irAf c H /~f 2 ~ 10 7 / c day -1 , 
for 7 = 10~ 8 . Now that we know the rate of events as 
a function of I, 9 and D C1 we can also evaluate the flux 
density in these terms 



s(i,e,D c ) = 



2ni 2 l 2 
D 2 (At) ohs 



e -u (l+z)w3 ( u (l + Z)W 3 > 1 

(12) 



and we assume a "top hat" cut-off at small angles, 
S(l,9,D c ) = 2ire- 1 ill 2 /D 2 c (At) ohs if w (l + z)l9 3 < 1. 
The event rate relevant to a survey is 

N(> S) = C [ -- D 2 dD sm9d9 (13) 

Here the integration volume V is constrained by the re- 
quirement that the flux be greater than S and that the 
observed duration of the event, (Ai) b s , be in the range 
in which the search was carried out. In our case, (A£) bs 
is in the interval (1 ms, 1 s) to coincide with the search 
parameters in Ref. [l[, and we will assume a standard 
flat cosmology. After integrating over 9 by hand, the 
remaining integrals were evaluated numerically, giving 



N(> 30 Jy) ~ 10°/c day" 



(14) 



whereas the single observed event gives a 99% double- 
sided confidence level estimate between 2 to 3 x 10 3 events 
per day [2j. (Based on the dynamics of smooth loops, one 
expects f c ~ 1.) The absence of fainter events, but still 
above the threshold of 0.3 Jy, is not significant since the 
prediction is only around 10 events above 0.3 Jy, and 
also the reach of the Parkes survey at threshold is less 
by a factor ~ 10 than at 30 Jy 0. This shows that the 
superconducting string model can give an event rate con- 
sistent with the Lorimer et al observations. Further, it 
is possible to show that for fixed observational frequency 
and duration range, N(> S) oc S~ x / 2 for large S. This 
is quite different from the S~ 3 / 2 fall off expected from 
uniformly distributed, identical sources. The slow fall off 
implies that there may be a population of very strong 
sparks e.g. ~ 1 per year above 1 MJy. 

A more thorough analysis, taking into account the con- 
straints imposed by different surveys, may impose further 
limits on the model. For example, the STARE survey [2(| 
places an upper bound of 7.5 x 10~ 2 events per day [2j for 
S > 80 kJy at an observational frequency of 61 1 MHz and 
durations ranging from 125 ms to a few minutes. With 
these observational parameters and the superconducting 
string model parameters used above, the model predic- 
tion is ~ 6 x 10 -2 events per day. This suggests that 
useful constraints on the string model may be placed by 
presently existing data (see Fig. [2]). In particular, a value 
of 7 significantly less than 10 -8 is already ruled out. 

While we expect a quasi-cusp to repeat every oscilla- 
tion, we do not expect the beaming direction to remain 
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FIG. 2: Event rates for PARKES and STARE survey param- 
eters. The vertical line at 30 Jy corresponds to the inferred 
rate from the PARKES survey and at 80000 Jy shows the 
constraint from the STARE survey. 



the same in every repetition. The angular momentum 
carried off by the beamed radiation is ~ Eb eam L and 
the moment of inertia of the string in the cusp region 
is ~ (ficr)L 2 , where a ~ (ujqL)^ 1 ^ 3 L is the length of 
string in the cusp region relevant to the observed radi- 
ation frequency. The angular velocity due to radiation 
backreaction is 9 ~ E\y eaim L/iiaL 2 . Since the time in- 
terval between two cusp events is ~ L, we estimate the 
change in the beaming angle as 



AO 



[iuL 2 



L 



V 



{lu^L) 1 ' 3 - 2 x 10~ 4 (15) 



which corresponds to a change in the observed flux 
AlnF| ~ \auj o L0l\ x 3A9/9 ~ 5 Hence the beamed 
fluence may get reduced by a factor e~ 5 ~ 10 ~ 2 . This 
estimate suggests that sparks due to consecutive cusps 
on strings would be beamed in only slightly different 
directions but, since the fluence is highly sensitive to 
the beaming direction (Eq. |T|)). the repeated event may 
or may not be observed, depending on the precise cusp 
shape and other parameters. In the 20 days that Lorimer 
et al searched for bursts, the particular loop would have 
on order 10 6 cusps, and since the solid angle of the beam 
is order 10 -6 , we can expect on order one event in those 
20 days. Also, cusps that occur much later, will occur 
with different parameters since the length and shape of 
the loop are likely to change due to backreaction. 

To model the observations, we have taken strings 
whose tension scale is rj ~ 10 14 GeV (so that 7 ~ 10~ 8 ) 
but with currents of only ~ 10 5 GeV. The small current 
compared to the tension scale can have a natural origin 
in terms of scattering of counter-propagating particles on 
the string 2l|,l22| and/or the modified dispersion relation 
due to background magnetic fields 23]. 

We have considered only the direct emission of 1.4 
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GHz radiation from the cusp. This corresponds to a 
very high harmonic emitted from the oscillating loop. 
The fundamental frequency is given by L^ 1 ~ 1 Hz, 
and 1 GHz emission corresponds to the 10 9 th harmonic. 
What happens to emission at lower frequencies? Since 
the cusp event is not expected to be in a host galaxy, 
it is surrounded only by the IGM with free electron 
number density n e ~ 10~ 7 cm~ 3 and plasma frequency 
Ljp = (47rn e e 2 /TO e ) 1 / 2 ~ 30 s" 1 or v p ~ 5Hz. The energy 
emitted from the cusp at v < 5 Hz cannot propagate 
in the IGM and must push it around, creating shocks, 
fireballs, and possibly gamma ray bursts (HI. Irl l24j. As- 
suming that all the very low frequency emission gets con- 
verted into gamma rays due to plasma effects, the power 
emitted is ~ %^ = 10 10 GeV 2 . For gamma rays at lGeV 
and for loops at a distance of lGpc, this gives a photon 
flux of 10~ 22 cm~ 2 s _1 which is far below the threshold, 
~ 10 _8 cm~ 2 s -1 of the Third Interplanetary Network. 

The electromagnetic emission from strings also distorts 
the cosmic microwave background spectrum and this ef- 
fect has been used to constrain superconducting strings 
2J, |25j. Our choice of parameters, G/i <C 10~ 6 , is within 
these constraints. It is worth pointing out that the earlier 
work on gamma ray bursts from cusp events primarily fo- 
cussed on superconducting strings placed in the galactic 
environment. More detailed investigation is needed to 
determine if gamma ray bursts are expected to accom- 
pany radio sparks occurring in the IGM. 

The superconducting string model may be tested in 
a variety of ways. Gravitational effects of strings with 
G/i ~ 10~ 10 will be weak and may not be within forsee- 
able detection capabilities, except possibly for gravita- 
tional wave bursts from cusps [26| . A promising possibil- 
ity is to look for signatures of particle emission, such as 
positrons [27I ] or other decaying particles produced where 
the current on the string quenches. These particles would 
give a distinctive feature in the emission from the vicin- 
ity of the spark. In addition to the observed sparker like 
events, there should be rarer events where we are even 
closer to the beam, such that ojL8 3 < 1. Then the spec- 
trum will not decay exponentially, and the characteristics 
of the event should be quite different. Also, kinks on su- 
perconducting strings will radiate in unusual "fan-like" 
patterns [28}. Most immediately, it is necessary to find 
more sparks and check if they are associated with host 
galaxies since supermassive strings of the kind we have 
considered are expected to roam outside of galaxies. 
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